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We present a dynamical mean-field theory (DMFT) study of the charge and orbital correlations 
in finite-size Lao.sCao.sMnOs (LCMO) nanoclusters. Upon nanostructuring LCMO to clusters of 
3 nm diameter, the size reduction induces an insulator-to-metal transition in the high-temperature 
paramagnetic phase. This is ascribed to the reduction in charge disproportionation between Mn 
sites with different nominal valence [Daa et al, Phys. Rev. Lett. 107, 197202 (2011)]. Here we 
show that upon further reducing the system size to a few-atom nanoclusters, quantum confinement 
effects come into play. These lead to the opposite effect: the nanocluster turns insulating again and 
the charge disproportionation between Mn sites, as well as the orbital polarization, are enhanced. 
Electron doping by means of external gate voltage on few-atom nanoclusters is found to trigger a 
site- and orbital-selective Mott transition. Our results suggest that LCMO nanoclusters could be 
employed for the realization of technological devices, exploiting the proximity to the Mott transition 
and its control by size and gate voltage. 

PACS numbers: 71.27.-l-a,71.10.Fd,71.30.+h,75.47.Gk 


I. INTRODUCTION 


Research on manganites dates back to the 1950’s, when 
Jonker and van Santeit^ reported the existence of fer¬ 
romagnetic metallic phase in mixed crystals of man¬ 
ganese oxides LaMnOa-CaMnOa, LaMnOa-SrMnOa, and 
LaMnOa-BaMnOa. However, the interest of a wide por¬ 
tion of the scientific community was only raised in the 
1990s, due to the experimental observation of a colossal 
magnetoresistance (CMR) effecti^i^ Indeed, the relative 
change in resistivity upon the application of an external 
magnetic field was much higher than the one observed in 
artificial magnetic/non-magnetic multilayer systems: up 
to 60% at room temperature in thin films.— Triggering 
such a CMR, however, requires cooling below the Curie 
temperature Tc and the application of relatively strong 
magnetic fields, preventing the technical application of 
the CMR effect to this day. Alternative routes to achieve 
CMR have also been followed in mixed valence mangan¬ 
ites such as the half-doped Lao.sCao.sMnOg (LCMO)^ 
which is insulating below 155 K and displays antiferro¬ 
magnetic and charge order often also accompanied by 
orbital order. Indeed, the antiferromagnetic insulating 
state is prone to instabilities. The transition toward a 
ferromagnetic metallic state can be triggered upon apply¬ 
ing a magnetic field^ doping, biaxial strain, pressure^ or 
an electric field— The experimentsiSril suggest that the 
destabilization of the charge-orbital order can also be ob¬ 
tained upon size reduction. In fact, it was observed also 
in Ndi_xCaa;Mn 03 ^d 2 and Smi_ 2 ;Caa;Mn 03 i^ in both 
at half-doping and in asymmetrically doped samples, as 


well as in Pro.sCao.sMnO^i^Ti^i compounds. Nonetheless, 
there remains some controversy^^ arising from the exper¬ 
imental difficulty to disentangle the effects of size reduc¬ 
tion from other effects, such as oxygen non-stoichiometry, 
disorder, strain etc. 

On the theoretical side, the effects of size reduction has 
been studied by density functional theory (DFT)-|-U22ri2i 
and DFT-|-dynamical mean-field theory (DMFT)— ^ The 
theoretical analysis is in remarkable agreement with the 
experiments and shows a correlation-driven destabiliza¬ 
tion of the charge-orbital order in bulk LCMO upon size 
reduction. This bears the prospects that, for the opti¬ 
mized size, a much smaller magnetic field is sufficient to 
trigger a CMR, as LCMO nanoclusters can be tuned to 
the verge of a metal-insulator transition. Hitherto, the 
DFT-I-DMFT calculations on this topic, e.g., those re¬ 
ported in Ref. were performed for bulk model systems 
with ab-initio parameters. That is, following the DFT 
calculations for 3 nm clusters, model bulk systems were 
constructed having the same unit cell volume as well as 
octahedral distortion as in the core of 3 nm cluster. This 
way, it was possible to take into account the interplay 
between strong electronic correlations within DMFT and 
the structural distortions induced by size reduction, ob¬ 
tained ab-initio through atomic relaxation within DFT. 
In this paper, we take a significant step forward, in terms 
of carrying out nanoscopic DMFT calculations for few- 
atom clusters with a DFT-derived tight-binding Hamil¬ 
tonian. This gives us the opportunity to consider the 
effect of size reduction from bulk to intermediate-sized 
clusters to few-atom clusters. Our calculations show an 
interesting evolution from the high-temperature param- 
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agnetic insulating (PI) state in the bulk to a paramag¬ 
netic metallic (PM) state in intermediate-sized clusters to 
re-entrant insulating solution for few-atom clusters. We 
also investigate the effects of applied hydrostatic pressure 
in the bulk, which turn out to be different than the ef¬ 
fects of size reduction. Considering few-atom clusters, we 
also show that electron doping, through the application 
of external gate voltage, drives an unexpected site- and 
orbital-selective Mott transition. 

The paper is organized as following. In Section |TT] we 
discuss the model employed for the description of mixed- 
valence manganites, and we describe the strategy we fol¬ 
lowed to include structural, finite-size, and many-body 
effects in the framework of a combined DFT-I-DMFT ap¬ 
proach. In Sec. m we present the DMFT results ob¬ 
tained for LCMO nanoclusters of different size. In par¬ 
ticular, in Sec. IIII Al we focus on the effects of quan¬ 
tum confinement on the spectral properties, while in Sec. 
IIII Bl we explore the effect of electrostatic doping by ap¬ 
plying an external gate voltage to the few-atom clusters. 
Finally, in Sec. IIVI we present our conclusions. 


II. METHOD: DFT+DMFT APPROACH FOR 
LA 0 . 5 CA 0 5 MNO 3 NANOCLUSTERS 

A. Bulk crystal and electronic structure 

Manganites, Ri-xAajMnOs with R being a trivalent 
rare-earth-metal element and A a divalent dopant, have 
a perovskite lattice structure, with the rare-earth atoms 
at corner positions, the Mn atoms at body center po¬ 
sitions and oxygen atoms at the face centered positions. 
Depending on the sizes of R and A, given by the so-called 
tolerance factor, the MnOe octahedra can tilt and rotate 
reducing the symmetry of the perovskite lattice from cu¬ 
bic to orthorhombic. With nominal oxygen valency 0^“, 
in half-doped compounds (i.e., x = 0.5) the manganite 
atoms are in a mixed valent Mn^-^+ state. This can lead 
to a charge disproportionation between the Mn sites in 
bulk half-doped manganites; in the extreme case one has 
a 50% of Mn"^+ sites with a configuration, and the 
other 50% of Mn^+ sites with a 3c?^ configuration. The 
charge-ordered state is associated with a real space or¬ 
dering of Mn^+/Mn^+ species in a 1:1 pattern. The oc¬ 
tahedral crystal held surrounding the Mn ions splits the 
3d orbitals into three low energy t 2 g orbitals {dxy, dxz, 
and dyz) and two higher energy tg orbitals (d 322_^2 and 
dx^-y^)- Due to the strong Hund’s exchange coupling, 
the three Mn electrons of Mn^’*' occupy the t 2 g orbital, 
in a high-spin St2g = 3/2 conhguration. The consequent 
energy gain associated to this state, makes it energeti¬ 
cally favorable for the extra electron of Mn^+ to occupy 
the doubly degenerate Cg states. Hence, the charge or¬ 
dering accommodates a static Jahn-Teller distortion at 
the Mn^+ sites, removing the degeneracy and lowering 
the symmetry of the system to monoclinic. 

The Mn Cg orbitals are delocalized due to a strong 


hybridization with the 02p states. On the contrary, Mn 
t 2 g orbitals do not hybridize strongly with the O and are 
localized. The magnetism is thus governed by the double 
exchange mechanism.— 


B. DFT calculations and downfolding procedure 

The first step of our study is a DFT calculation of 
both LCMO bulk and nanoclusters. To this end, we 
use projected augmented wave (PAW) pseudopotentials 
with an energy cutoff of 450 eV and performed calcu¬ 
lations within a spin-polarized generalized gradient ap¬ 
proximation (GGA)^ as implemented in the Vienna ab- 
initio Simulation Package (VASP)— The forces on the 
atoms are converged to less than 0.01 eV /A. All DFT cal¬ 
culations for bulk LCMO was performed with a 4 x 2 x 4 
/c-mesh. On the other hand, for the nanoclusters we 
use only the F point for the /c-space integration. The 
structural optimization of the bulk LCMO is performed 
considering both antiferromagnetic (AFM) CE and fer¬ 
romagnetic (FM) configurations. We find the AFM CE 
phase to be lower in energy compared to the FM phase 
by 45 meV/f.u. The stoichiometry of the system is main¬ 
tained in the construction of the nano-clusters, as dis¬ 
cussed in Refs. [H and [H. 

A necessary step for the subsequent DMFT calculation 
is the extraction of the relevant tight-binding parameters, 
calculated ab initio within DFT. To achieve this, we em¬ 
ploy the downfolding method as implemented in the N-th 
order muffin tin orbital (NMTO) basis^i^ with poten¬ 
tial parameters borrowed from self-consistent linearized 
mufRn-tin orbital (LMTO)^ calculations. Through the 
NMTO-downfolding procedure, a low-energy Hamilto¬ 
nian H{k) involving only Mn Cg Wannier orbitals^ is 
constructed in fc-space by integrating out all other de¬ 
grees of freedom. The Fourier transformation of H(k) 
provides the tight-binding parameters. In the follow¬ 
ing we refer to this structure as S'buik- We use S'buik to 
describe bulk manganites within standard DFT-(-DMFT 
calculations. In order to construct a low-energy Hamil¬ 
tonian for the nanoclusters, as explained in Ref. [ 2 ^ a 
cluster of nearly spherical shape having 3 nm diameter is 
cut out from a large supercell of the bulk crystal structure 
in monoclinic P2i/m symmetry, which is then subject to 
a full structural optimization. In the following we refer to 
this structure as ^nano- The NMTO-downfolding calcula¬ 
tion is then carried out on the self-consistent LMTO cal¬ 
culation for a model bulk structure, referred to as S'modei, 
which is constructed considering the MnOe octahedra se¬ 
lected from the core region of S'nano, and applying various 
symmetry operations, as explained in detail in Ref. [ 2 ^. 
The DF-derived tight-binding Hamiltonians contain the 
information of the structural as well as electronic changes 
at the level of one-electron theory that happen upon size 
reduction. 

In order to take into account many-body effects within 
DMFT, we build a low-energy model (discussed in detail 
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in Sec. Ill Cl) with the DFT-derived hopping parameters 
as an input. However, there is fundamental difference 
between the calculations performed in this work and that 
of Ref. m as explained is the following: In Ref. we 
performed DMFT calculations on the constructed model 
bulk structure (5'inodei)- Instead, in the present work, we 
use the same parameters of S’modei to construct few/atom 
nanoclusters of different size, which we solve within the 
nano-DMFT scheme discussed in Sec. mm 


C. Low-energy effective eg + St 2 g model 


The low-energy Hamiltonian describing the mangan- 
ites is given as below 


ijmm' aa' im 


+u E 




im<7 


imm' (j(j 


( 1 ) 


Here, denotes the generic matrix element of the 

one-particle DFT Hamiltonian in the basis of the the 
downfolded NMTO Wannier orbitals m, m' at site i, 
j] the on-site Coulomb interaction between the Cg elec¬ 
trons is parametrized in terms of an intra-orbital repul¬ 
sion U = 5 eV, a Hund’s exchange J = 0.75 eV, and an 
inter-orbital interaction U' = U — 2J The values of 
the interaction parameters, taken from the literature^ 
represent realistic estimates for manganites. Further¬ 
more, the Cg orbitals are coupled to a (classical) disor¬ 
dered spin S' by a Hund’s exchange J. This term rep¬ 
resents the interaction between the itinerant Cg electrons 
with the localized electrons in the half-filled t 2 g mani¬ 
fold. In the notation adopted here, S = ±1 is the classical 
spin corresponding to the high-spin states of St^g =±3/2 
while its modulus is included in the value J7 = 1.35 eVi^ 
Hamiltonian ([1]) represents a standard low-energy model 
for the description of electronic correlations in mangan- 
ite compounds^^— Recently, the classical spin descrip¬ 
tion of the t 2 g manifold has been thoroughly revisited4S 
Thanks to the significant technical improvements in the 
field of impurity solvers, this allows the direct treatment 
of a hve-orbital model for the whole Mn3c? multiplet, in¬ 
cluding also spin-flip and pair-hopping terms beyond the 
density-density Coulomb interactionj ^^d^ Remarkably, in 
the case of pure LaMnOa, a detailed analysis showed 
an excellent agreement between a classical spin and a 
full quantum many-body treatment of the t 2 g orbitals, 
e.g., for the Cg spectral functions4£ Moreover, the es¬ 
timate obtained for the spin-spin correlations functions 
{St” Sl^”) ^0.74 was found to be consistent with the pic¬ 
ture of aligned Cg and S = 3/2t2g spins4S In light of these 
considerations, the restriction to the low-energy model 
described by Hamiltonian m represents a realistic and 
physically sensible choice to study correlation effects in 


Mn compounds with a half-filled t 2 g manifold and par¬ 
tially filled Cg orbitals. 


D. Dynamical mean-field theory with inequivalent 
Mn atoms in the unit cell 

In the following we discuss the technical details for the 
solution of the many-body Hamiltonian o in the frame¬ 
work of DMFT— for inhomogeneous systems— “— 
That is, we solve an auxiliary Anderson impurity prob¬ 
lem for each inequivalent Mn site in the unit cell. More¬ 
over, we perform an average over the disordered classical 
t 2 g spin S'— This procedure yields a local 2x2 self¬ 
energy in the Cg manifold of each of the Mn atoms in the 
unit cell, while neglecting non-local self-energy elements 
between different Mn atoms. Let us start by discussing 
the bulk DFT-I-DMFT calculations, which have been per¬ 
formed on Sbuik and Smodei structures to describe LCMO 
bulk and nanoscopic (3 nm) clusters. The bulk mono¬ 
clinic unit cell of LCMO contains eight Mn atoms but 
only three kinds of Mn atoms are locally inequivalent, 
labeled as Mnl(l), Mnl(2), and Mn2. The four Mnl 
atoms have a nominal valence 3-1- and occupy the bridge 
sites of the zig-zag ferromagnetic chain in the CE-type 
magnetic order that sets in below the Neel temperature 
in bulk LCMO. Those Mnl are further divided into two 
Mnl(l) and two Mnl(2) sites, by symmetry. The four 
Mn2 atoms have a nominal valence 4-|- and occupy the 
corner sites of the zig-zag ferromagnetic chain. Exploit¬ 
ing these symmetries, the overall computational effort for 
the bulk amounts to the solution of three auxiliary im¬ 
purity problems within a DMFT self-consistent scheme, 
with the corresponding subtraction of three inequivalent 
double counting terms in 

= ( 2 ) 

where U = U — denotes an averaged interaction and 
are the NMTO orbital occupancies for each of the 
inequivalent Mn in the unit cell.— In the numerical 
calculations, we employ a Hirsch-Fye Quantum Monte 
Carlo^ impurity solver, with a Trotter discretization 
Ar^ Ri 0.027 and inverse temperature /3 = 20 eV“^ in 
the paramagnetic phase. 

For the nano-DMFT calculations, we consider the 
LCMO finite-size nanoclusters shown in Fig. [T] and de¬ 
scribed by the structure 5nano- The symmetry of the 
nanoclusters is much lower than in the bulk due to finite- 
size effects. This leads to many more inequivalent Mn. 
For example, the TV = 46 Mn atom cluster contains 
TVineq = 23 inequivalcnt Mn sites, as it possesses only 
the inversion symmetry with respect to the center of the 
cluster. In this case, we need to solve for TVjneq impu¬ 
rity problems, yielding TVjneq 2x2 local self-energy ma¬ 
trices From these, new (TV x 2) x (TV x 2) cluster 

Green’s functions are calculated by solving the Dyson 
equation, enforcing the self-consistency at the level of 
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FIG. 1: (Color online) Schematic representation of the Mn 
sublattice of the few atom sized small clusters considered in 
the present DMFT work, having N = 46, 20, 8, and 4 Mn sites. 
The clusters are built by chopping off atoms from the super¬ 
cell described by the tight-binding Hamiltonian of Snano. The 
red (gray), blue (dark gray) and green (light gray) atoms cor¬ 
respond to Mnl(l), Mnl(2), and Mn2, respectively, according 
to their classification in the bulk. The cluster boundaries cor¬ 
respond to dangling Mn-Mn bonds. La, Ca, and O atoms 
(not shown) are effectively taken into account in the ab-initio 
parameters of 5'nano through the downfolding procedure. 


the whole nanocluster. This scheme has been employed 
for all the few-atom Mn clusters, shown in Fig. [TJ taking 
into account the specific symmetries of each structure. 
We stress once more that this nano-DMFT approach is 
different in spirit from that adopted in Ref. l22l in which 
bulk-DMFT calculations were performed, and the effect 
of size reduction was considered via the change in tight- 
binding parameters of S'modei with respect to S'buik- In¬ 
stead, within nano-DMFT calculations we explicitly take 
into account the boundary effects of the finite-size nan¬ 
oclusters constructed with the tight-binding parameters 
extracted from 5'nano- 


III. RESULTS FOR FINITE SIZE LCMO 
NANOCLUSTERS 

A. Spectral properties, charge and orbital order 

In the following we discuss the spectral properties of 
the nano clusters shown in Fig. 1. The DMFT spectral 
functions are shown in Figs. [2J for three representative 
Mn atoms belonging to Mnl(l), Mnl(2) and Mn2 kinds, 
following the bulk classification (even though for the nan¬ 
oclusters the atoms of e.g. Mnl(l) type are of course not 
equivalent any longer as in 5modei)- At the outset, we 
notice that the size reduction has profound consequences 
on the electronic structure of the LCMO nanoclusters, 
in the sense that the metallic character decreases grad¬ 
ually in moving from = 46 to 4. In the extreme case 
of A^ = 4, the system is insulating even at the high tem¬ 
perature considered here (/3 = 20 eV“^), reflecting the 
strong quantum confinement effects induced upon size 
reduction. The observed metal-to-insulator transition is 
accompanied by a strong enhancement of the charge dis¬ 


proportionation between Mn'^+ and Mn^+. We further 
find that size effects are accompanied by an overall en¬ 
hancement of the orbital polarization. These effects are 
strong for Mnl(I) and Mn2 atoms which are located at 
the core of the clusters. Instead, all Mnl(2) sites are 
located at the surface of the clusters (with maximum 
number of dangling bonds) and have a nearly insulating 
and almost fully orbitally-polarized spectral function ir¬ 
respective of cluster size. 

In the bulk LMCO, the orbital polarization between two 
Cg orbitals, 3z^—and x^—y^ is complete for Mnl atoms 
with nominal 3-|- valence, and zero for Mn2 atoms with 
nominal 4-|- valence. Moving to cluster of 3 nm diame¬ 
ter it was shown^^ that both charge disproportionation 
between Mnl and Mn2 atoms, as well as the orbital po¬ 
larization at the Mnl atoms decreases considerably com¬ 
pared to the bulk, leading to metallicity. Indeed, we find 
that upon further reduction of the size to few-atom nan¬ 
oclusters, quantum confinement effect comes into play, 
making the charge disproportionation and orbital polar¬ 
ization increase again, especially for Mnl(l) atoms. Both 
charge disproportionation and orbital polarization show 
an increasing trend upon decreasing the cluster size, from 
= 46 to 20 to 8 until it becomes almost complete for 
the TV = 4 nanocluster. 

In order to study quantitatively the effects of size re¬ 
duction, we consider the cluster averaged charge order 
and the orbital polarization. To this end, we define the 
occupation rii and the polarization pi of Mn site i as 

15} ^ 3 ^ 

Pi = 2 

{S} rna 


where the average over {5} takes into account the two 
possible configurations of the classical t 2 g spin 5 (as dis¬ 
used in Sec. IIIDI) . Hence, using the definitions above, we 
compute the following quantities 


(An) 

(Apa) 


1 

T^ 


5]n, - 

i^cx. 


1 


i^cx. 


1 


ie/3 


( 4 ) 


where a ^ fi denotes Mnl and Mn2 sites, respectively, 
while TVq,, Np denote the number of Mn sites of the cor¬ 
responding kind in the nanocluster. In particular, (An) 
is the charge disproportionation averaged over the nan¬ 
ocluster, while (Apa) is the orbital polarization averaged 
over all Mn site in the nanocluster belonging to the same 
kind a. Note that in Eq. the occupation n^ and the 
polarization pi for the different Mn sites i (even those 
belonging to the same Mn kind a) are in general inequiv¬ 
alent because of the symmetry of the finite nanoclusters 
constructed from 5nano, which are different from the bulk 
models described by 5buik or 5niodei- 
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FIG. 2: (Color online) Evolution of the spectral properties of the nanoclusters, shown in Fig. 1, described by the tight-binding 
parameters of Snano. Each row shows the spectral function A{u:) for a cluster having N Mn atoms. The red (gray), blue (dark 
gray) and green (light gray) lines correspond to representative Mnl(l), Mnl(2), and Mn2 atoms, respectively, indicated by 
shaded spheres in the nanoclusters shown in Fig. [T] The dark and light filled curves in each panel denote the contributions 
from the and orbitals, while the black solid line denotes the on-site spectral density of the Cg manifold. A 

metal-to-insulator transition (MIT) upon decreasing cluster size, due to quantum confinement effect, is evidenced accompanied 
by an enhancement of the charge and orbital order. 


The cluster-averaged (An) and (Apa) are shown in 
the right-hand side of each panels in Fig. [^as a function 
of the (inverse) number of the Mn atoms in the cluster. 
For comparison, on the left-hand side of each panel the 
results obtained in Ref. [12 within DFT-f-DMFT calcula¬ 
tion performed on the bulk nanomodel Smodei are also 
shown (applicable to 3 nm diameter cluster containing 
about 200 Mn atoms). As a general trend, the charge 
order and orbital polarization are found to be the largest 
for the smallest cluster and decrease upon increasing the 
cluster size. In particular, An shows a smooth reduc¬ 
tion with increasing system size from a value close to 
An = 1 for A = 4, corresponding to the limit in which 
the Mn2 sites are completely empty, to An « 0.70 for 
N = 46. Those values can be compared to An ~ 0.35 
found for the 5'nano model. A similar behavior is found 
for the averaged orbital polarization, although the be¬ 
havior is somewhat dependent on the system shape and 
symmetry, especially for the smaller nanoclusters consid¬ 


ered here. However, in general the orbital polarization 
of the Mnl atoms tends to decrease, by increasing clus¬ 
ter size, while that for the (almost empty) Mn2 atoms is 
always negligible. A more careful analysis (not shown) 
reveals that the averaged Ap is systematically larger than 
the one of the core Mn atoms (especially for larger clus¬ 
ters), as it also includes the contribution of Mn atoms 
on the surface of the cluster, which are characterized by 
sharper spectral structures and a stronger orbital polar¬ 
ization. 

If we also consider the results obtained from the 
DFT-I-DMFT calculations of Ref. [12, we find that both 
the charge order (An) and the orbital polarization (Apa) 
are strongly enhanced for 5buik with respect to the 3 nm 
cluster size described by iSmodei, which drives a metal-to- 
insulator transitioni^ 

The above observations are depicted by the following sce¬ 
nario: Starting from the smallest size cluster and upon 
progressively increasing the system size, one encounters 
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FIG. 3: (Color online) Cluster averaged charge order (An) 
(upper panel) and orbital polarization (Apa) (lower panel) 
of LCMO bulk and clusters. On the right-hand side of each 
panes we show the nano-DMFT results obtained for the finite- 
size clusters built with the parameters of Snano, as a function 
of the (inverse) number of the Mn atoms in the cluster. On the 
left-hand side of each panel we also show the DFT-I-DMFT 
results obtained within the calculations of Ref. for the bulk 
model Smodei (applicable for 3 nm diameter cluster containing 
about 200 Mn atoms) and on bulk LCMO S'buik. Those are 
compared to additional DFT-I-DMFT result for Spress, cor¬ 
responding to the bulk under the application of hydrostatic 
pressure, resulting into 6% volume reduction. The light brown 
(light tray) shade in the upper panel represents the transitions 
between PI and PM phases. Arrows and dashed lines are a 
guide to the eye. 


TABLE I: Orbital DFT-I-DMFT occupancies for the three in¬ 
equivalent Mn atoms in the unit cell of Smodei and Spress. 
In brackets we give the corresponding occupancies for the 
one-particle low-energy DFT Hamiltonian without the effect 
of DMFT correlations, which strongly enhance the charge- 
orbital order. Both the orbital polarization in Mnl atoms, as 
well as the charge disproportionation between Mnl and Mn2 
atoms are more pronounced for the nano model Smodei than 
for the pressurized bulk system Spress. 



Pressure 6% 

Nano model 


3«^ — 

3«^ — 

Mnl(l) 

0.52 (0.39) 0.06 (0.11) 

0.52 (0.31) 0.09 (0.20) 

Mnl(2) 

0.36 (0.34) 0.07 (0.13) 

0.72 (0.38) 0.04 (0.19) 

Mn2 

0.33 (0.31) 0.17 (0.21) 

0.16 (0.21) 0.16 (0.25) 


an insulator-to-metal transition in few-atoms nanoclus¬ 
ters, which is driven by a weakening of quantum con¬ 
finement effects and charge-orbital correlations. With 
increasing system size, the nanocluster smoothly evolves 
towards the results obtained with the bulk nanomodel, 
which is metallic due to the weak structural distortions 
included in the parameters of Smodei- On the other hand, 


bulk LCMO Sbuik is strongly distorted. The structural 
distortions lead to the enhancement of charge-orbital 
order and drive the system across a second metal-to- 
insulator transition between Smodei (3 nm) and Sbuik 
(bulk),^ 

It is interesting to compare the results obtained for 
Smodei and for the bulk under applied hydrostatic pres¬ 
sure Spress corresponding to the same volume reduction 
of the unit cell, i.e., 6% for 3 nm nanocluster described by 
Smodei- We find that the effects of hydrostatic pressure 
and the size reduction are rather different. The results for 
Spress indicate that charge order and orbital polarization 
are strongly reduced with respect to the case of Smodei , as 
shown in detail in Table HI In particular, we find that for 
Spress the charge disproportionation between Mn^+ and 
Mn^T is negligible, and it is accompanied by an overall 
reduction of the orbital polarization. This makes the sys¬ 
tem under hydrostatic pressure to be far more metallic 
compared to a 3 nm cluster. 


B. Site- and orbital-selective Mott transition 
driven by applied gate voltage 


The interesting insulator-to-metal-to-insulator transi¬ 
tion in LCMO upon size reduction is based on a complex 
interplay between quantum confinement effects and the 
structural distortions occurring upon size reduction from 
the bulk LCMO to 3 nm nanoclusters. As the onset of 
the peculiar charge- and orbital-ordered state found in 
the bulk LCMO relies on the balance between Mn^+ and 
Mn'^+j it is interesting to study the effect of electron dop¬ 
ing on the few-atom clusters. A change in the number 
of carriers without changing the chemical composition of 
the system can be achieved by the application of an ex¬ 
ternal gate voltage Vg. In this section, we investigate this 
issue through DMFT calculations considering the clusters 
with TV = 46 and 4 Mn atoms. We consider the limiting 
case where there is only an inhnitesimally small tunneling 
contact with the environment, so that the we can account 
for the gate voltage by changing the DMFT chemical po¬ 
tential. We neglect the effect of doping on the DFT effec¬ 
tive potential, as in the virtual crystal approximation!^ 
In Fig. |4] we show t he orbi tally-resolved spectral weight 
at the Fermi energy A{Ef) (averaged over a unite energy 
window ^ T) as calculated from the Green function at 
imaginary time r = ,5/2, and the occupations of the two 
Cg orbitals for representative sites of ^nano with TV = 46 
and TV = 4 Mn sites clusters. The value Vg = 0 corre¬ 
sponds to the results in the previous sections and to an 
average cluster occupation (n) = 0.5 electrons in the Cg 
orbitals. Upon changing Vg, we increase the number elec¬ 
trons in the LCMO cluster (electron doping). When the 
low-lying Cg orbital (e.g., the x"^ — in the case of Mnl) 
on a Mn site becomes half-filled, strong electronic corre¬ 
lations drive an orbital-selective metal-to-insulator tran¬ 
sition, with the opening of a Mott gap. Such an orbital- 
selective Mott transition has been reported before for the 
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FIG. 4: (Color online) Evolution of the spectral weight at the Fermi energy A{Ef) (averaged over a unite energy window ~ T) 
and of the orbitally-resolved occupation with gate voltage Vg (electron doping) for the representative Mn sites of the = 46 
(upper panels) and = 4 (lower panels) clusters. The red (gray) blue (dark gray) and green (light gray) curves correspond 
to representative Mnl(l), Mnl(2), and Mn2 atoms, respectively, while the light and dark colors denote 2>z^ — and 
states. The gate voltage Vg drives an orbital selective Mott transition, which becomes sharper upon decreasing the system size. 


Hubbard model, originating from different bandwidth (or 
correlations) for different orbitals,— or due to the band 
degeneracy lifting,— as well as for materials.— The 
charge disproportionation between Mn^+ and Mn^+ also 
results in a strong site-selective character of the transi¬ 
tion. Site-selective behavior of similar kind has also been 
reported recently for bulk system s^^i^^ . In our case, we 
have one insulating orbital which is integer filled (n = I), 
but neither the occupation of the other (metallic) or¬ 
bital, nor the cluster average electron density is integer. 
One possible interpretation of this novel orbital- and site- 
selective Mott transition is associated to the role of the 
Hund’s exchange coupling to the Og-spins. A filling of 
one electron in an orbital thus means that all states with 
spin parallel to the t 2 g-spins are occupied, while those 
with opposite spin are empty. 

The orbital-selective MIT is relatively homogeneous in 
the cluster, meaning that Cg Mn sites of the same kind 
turn insulating at a similar value of Vg. On the other 
hand, the critical value of Vg still depends both on the 
kind of Mn site and the cluster size. Further increasing Vg 
increases the population of the other eg orbital (e.g., the 
3 ^ 2 _j ,2 of Mnl). Above Vg «2.6 eV (until the 

chemical potential lies within the Mott gap) all orbitals 


are half-filled and Mott insulating, with A{Ef) = 0 (see 
upper panels in Fig.|4|) and display a divergent imaginary 
part of the DMFT self-energy (not shown). In general, 
we observe that, upon decreasing the cluster size, the 
orbital-selective MIT is found at a smaller value of Vg. 
The transition also appears to be sharper upon changing 
Vg. This can be understood by considering the more lo¬ 
calized nature of the Cg orbitals and the enhanced orbital 
polarization observed for smaller cluster sizes. This ef¬ 
fect is important in view of possible applications. For an 
appropriate system size, half-doped LCMO manganites 
nanocluster can be driven across a MIT by applying an 
external gate voltage. 


IV. CONCLUSIONS 

We investigated the effects of size reduction on the 
charge and orbital order in the LCMO mixed valence 
manganite within the DFT-fDMFT framework. As was 
shown before,— the size reduction from bulk to a 3 nm 
nanoclusters weaken the distortions the bulk crystal 
structure and induces an insulator(bulk)-to-metal(3 nm) 



































transition in the high-temperature paramagnetic phase, 
along with a weakening of charge and orbital dispropor¬ 
tionation. Here, we extend the analysis by considering 
nanoclusters of just a few-atoms. Upon reducing the sys¬ 
tem size we observe the opposite trend: driven by the 
quantum confinement, there is a second metal(3 nm)-to- 
insulator(few atoms) transition and an enhancement of 
both charge and orbital disproportionation. We also in¬ 
vestigated the effect of electron doping on the few-atom 
nanoclusters by applying an external gate voltage. We 
observe an orbital-selective Mott transition at a critical 
value of the gate voltage which corresponds to an inte¬ 
ger filling of only an individual Cg orbital. The orbital- 
selective nature of the transition is a direct consequence 
of the orbital polarization and the strong Hund’s ex¬ 
change splitting. At the same time, the strong charge dis¬ 
proportionation between Mn^+ and Mn'^^ sites induces 
also a site-selective character, with different Mn kinds 
turning insulating at different values of the gate voltage. 

Our theoretical prediction of a reentrant insulator-to- 
metal-to-insulator transition and the reported gate volt¬ 
age control calls for further experiments. Technical ap¬ 


plications are discernible since these two control param¬ 
eters should allow us to fine tune LCMO nanoclusters to 
the verge of a Mott transition. In this situation smallest 
changes in temperature, voltage, magnetic field etc. can 
trigger a gigantic change in conductance. 
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